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Natural Hybrids between Ligularia vellerea and 
L. subspicata ( Asteraceae: Senecioneae ) 
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Abstract; Three types of intermediate morphological individuals were found in the sympatric distribution areas of 
Ligularia vellerea and L. subspicata. These intermediate individuals were hypothesized to be the hybrids of these two 
species based on the detailed comparison of their morphological characters. To test this hypothesis, we compared the 
DNA sequence of L. vellerea, L. subspicata and the intermediate individuals by using the internal transcribed spacer 
region of nuclear ribosomal and atpB-rbcL intergenic spacer region of the chloroplast. Direct sequencing of ITS clear- 
ly had additives and then cloned sequencing cloned out two kinds of sequences that were identical to those putative 
parental species, thus these putative hybrids were approved. With different morphological characters among three 
types, their plastid donors were all L. subspicata according to the maternal sequences of atpB-rbcL. 


Key words; Ligularia vellerea; L. subspicata; Hybrid progeny; ITS; atpB-rbcL 


There are increasingly voices to illustrate the 
important role of natural hybridization in plant evolu- 
tion (Arnold, 1997; Rieseberg, 1997; Rieseberg, et 
al., 2000; Abbott et al., 2008, 2010; Soltis and 





Soltis, 2009; Chase et al., 2010). The consequences 
of hybridization may include introgression affecting 
one or both taxa, formation of hybrid (especially al- 


loploid) species, and development of reticulate pat- 
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terns of evolution within a group ( Arnold, 1997; 
Soltis and Soltis, 2009). 

The genus Ligularia Cass. consists of about 140 
species mainly distributed in eastern Asia, one-hun- 
dred twenty-four species of which are known in Chi- 
na (Liu, 1989; Liu et al., 1994; Liu and Illariono- 
va, 2011). The majority of Ligularia distributed in 
the Hengduan Mountains area are endemic ( Liu et 
al., 1994; Liu et al., 2006). Interspecific diploid 
hybridization was considered to be an important 
cause for the high species diversity in this genus 
(Liu et al., 2001, 2006). Natural hybrids were 
commonly found in certain area which Ligularia 
sympatric distributed (Liu et al., 2006; Yu et al., 
2011). LigulariaXmaoniushanensis X. Gong & Y. Z. 
Pan is the first reported natural hybrid species 
through comprehensive studying ( Pan et al., 2008). 
Furthermore, by comparing L. lamarum and L. sub- 
spicata chemical similarities, Saito et al. (2011) 
presumed that the two species hybridized with other 
species of this genus commonly. Yu et al. (2011) 
studied one mix population with 6 Ligularia species, 
based on morphological and molecular data (cpDNA 
and ITS sequences, ISSR markers), bidirectional but 
asymmetric hybridization and also backcross was de- 
tected between L. subspicata and L. nelumbifonia. 
Hybridization is very common in Ligularia with vari- 
ous form and progenies. 

During field collection trips to the Hengduan 
Mountains in Yunnan, we found a mixed population 
including L. subspicata (Bur. et Franch.) Hand.-Mazz. , 
L. vellerea ( Franch.) Hand-Mazz., and a putative 
natural hybrid occurring on Mt. Daxue. The putative 
nothospecies was considered to be a natural hybrid 
between L. subspicata and L. vellerea based on its in- 
termediate morphological appearance and the partial 
overlap of the flowering period. Given this anthesis 
overlap, cross-pollination between two putative par- 
ents could occur naturally under sympatric conditions. 

In the current study, the ITS region of the nu- 
clear ribosomal DNA of all the sampled individuals 


was sequenced to determine the hybrid status of the 


morphologically taxon. Once its hybrid status was 
confirmed, we sequenced chloroplast atpB- rbcL re- 
gion for three taxa to determine the relationship be- 


tween these hybrids and the two natural species. 


1 Materials and methods 
1.1 Material collection and molecular methods 

We collected Ligularia vellerea, L. subspicata and 
the putative hybrids in the mix population on the Dax- 
ue Mountain, Shangri County, Yunnan, China. Young 
leaves of 25 individuals (including 10 L. vellerea, 10 
L. subspicata and 5 putative hybrids) were collected 
and dried in silica gel in field for DNA extraction. 
Though no other Ligularia species was found in this 
population, a set of other species (sequences down- 
load from NCBI, Appendix) were included to verify 
if this population genetic material contaminate by 
other species in the phylogenetic analysis. Total ge- 
nomic DNA was extracted from the silica gel-dried 
leaf tissue using the cetyltrimethyl ammonium bro- 
mide (Cruz-Mazo et al., 2009) methods (Doyle and 
Doyle, 1987) with minor modification. 

The ITS region of the nuclear ribosomal DNA of 
all the sampled individuals was amplified using 
primers ITS4 and ITS5 ( White et al., 1990). PCR 
was conducted in a total reaction volume of 20 wL, 
containing 10-30 ng template DNA, 2.0 uL 10x 
PCR Buffer with (NH,),50,, 1.3 uL MgCl, (25 
mmol +-L') , 0.5 uL dNTP (2.5 mmol +L" each), 
0.3 mmol + L™' each primer and 0.75 unit of Taq pol- 
ymerase (Takara). Amplifications were performed as 
following conditions: 1 cycle, 94 C , 5 min; 30 cy- 
cles, 94 C , 45 sec; 55 C, 1 min; 72 C, 40 sec, 
followed a final extension of 1 cycle, 72 °C , 5 min. 
The PCR products were purified by electro-phoresis 
through a 1.2% agarose gel followed by use of an 
E. Z. N. A. Gel Extraction Kit (Omega, Guangzhou , 
China). All accessions were subjected to sequen- 
cing, with amplification primers in an ABI 3700 
DNA automated sequencer with the BigDye Termina- 
tor Cycle Sequencing Kit ( Applied Biosystems , Fos- 
ter City, California, U.S. A.). The sequences were 
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aligned and compared in SeqMan ( DNASTAR, Bei- 
jing, China) (Burland, 1999). While direct se- 
quencing is feasible for L. vellerea and L. subspicata , 
it produced chimeric or unreadable peaks in the 
chromatogram for the putative hybrids. Hence, clo- 
ning sequencing was performed subsequently for the 
putative hybrid. Purified PCR products were cloned 
into plasmids using the pUM-T vector system 
( Bioteke). Four to twenty positive clones were se- 
lected for each amplification product and cultured for 
isolating plasmids. Positive clones with the inserts of 
correct size were confirmed by colony PCR. The 
plasmids with correct inserts were sequenced using 
ITS4 and ITSS primers. 

All samples from L. vellerea, L. subspicata and 
the putative hybrid were used in the chloroplast 
atpB- rbcL. Chloroplast atpB- rbcL region were am- 
plified using primers atpB and rbcL (Chiang et al., 
1998). PCR was conducted in a reaction volume of 
20 uL, containing 10-30 ng template DNA, 2.0 uL 
10xPCR Buffer with (NH,),50,, 1.0 uL MgCl, 
(25 mmol + L`), 1.0 pL dNTP (2.5 mmol - L'! 
each), 0.3 mmol + L”' each primer and 1.5 unit of 
Taq polymerase (Takara). Amplification conditions 
were as follows: 1 cycle, 94 C , 4 min; 30 cycles, 
94 C , 45 sec; 53 C, 45 sec; 72 C, 1 min or 50 
sec, followed by 1 cycle, 72 C , 7 min. PCR prod- 
ucts were purified and then directly sequenced bidi- 
rectional using the methods mentioned above. 

1.2 Phylogenetic analysis 

Minor variation (usually one or two nucleotide 
sites) between some clones may have resulted from 
PCR error caused by Taq DNA polymerase. Putative 
PCR-mediated recombinants were excluded before 
further phylogenetic analysis. All sequences were a- 
ligned in Clustal X (Thompson et al., 1997) com- 
plimented PAUP4.0. Information about variable 
sites was obtained using the program DnaSP 4. 0 
( Rozas et al., 2003). The aligned sequences were 
used to infer the phylogeny of the sampled individu- 
als using the criterion of maximum parsimony, im- 


plemented in PAUP * version 4. Ob (Swofford, 2002). 


Parsimony analyses were carried out using a heuris- 
tic search with tree bisection and reconnection 
branch swapping, the Multrees option, ACCTRAN 
optimisation, and 1 000 random addition replicates 


for the nuclear and chloroplast datasets. 


2 Results 
2.1 Morphology 

Ligularia vellerea and L. subspicata belong to 
Ser. Scapicaulis S. W. Liu and Ser. Ligularia Cass. , 
respectively (Liu, 1989). The main morphological 
differences between the two species are leaf shapes, 
hairs on stems and capitula types (Table 1, Fig. 1). 
The leaves of L. vellerea are ovate with cuneate base, 
whereas those of L. subspicata are cordate or hastate. 
The stem of L. vellerea is full of densely long white 
hairs, whereas that of L. subspicata is proximally gla- 
brous. There are 5—7 ray florets in one capitulum of 
L. vellerea , but none in L. subspicata. The putative nat- 
ural hybrids are intermediate between those of the 
putative parental species, with ovate or cordate leav- 
es, with/without ray florets in capitula and with/ 
without hairs on stem and the morphological charac- 
ters show three types of apparent transition state. 
Type-a and type-b are similar to L. vellerea, but 
there are thinner hairs on stem in type-a and type-b, 
and no ray florets in type-b. Type-¢ is similar to 
L. subspicata but with ray florets. The shape form of 
their leaves show apparently transition state. 
2.2 Phylogenetic analysis based on nrDNA 
ITS4-5 

The direct sequence aligned matrix generated a 
total of 642 bp or 641 bp (because of one indel). No 
variation was detected within species of L. vellerea 
and L. subspicata. There were 15 variation sites be- 
tween the two putative parental species. The putative 
natural hybrids had more than ten double peaks, and 
cloned sequences obtained two sequence types that 
were identical to those of L. vellerea and L. subspicata , 
respectively, which then clearly displaying that all 
double copies of ITS came from the putative parental 


species (Table 2). 
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L.vellerea H-a H-b H-c L. subspicata 


Fig. 1 From up to bottom show stems hairs, florets, laminas and whole individual of putative hybrids and their putative parents respectively. 
L. vellerea show the ovate with cuneate base leaves, 5-7 ray florets in one capitulum and full of densely long white hairs on stems; 
H-a show the cordate leaf, with stem leaves, several ray florets in one capitulum and thinner white hairs on stems; H-b show the 
cordate leaf, with stem leaves, no ray florets in capitulum and densely long white hairs on stems; H-e show the cordate leaf, 
with stem leaves, several ray florets in one capitulum and glabre stems; and L. subspicata show the hastate or cordate leaves 


with stem leaves, no ray florets in one capitulum and glabre stems. Photographs not to scale 


Table 1 Main morphological differences characterizing of Ligularia vellerea, L. subspicata and three types of putative hybrids 








Characters 
Taxa 
stems clothing stem leaves sexuality of florets composite raceme 
Ligularia vellerea densely long white puberulent without stem leaf numerous tubular florets, with several ray florets 
putative hybrid a ( H-a) sparsely long white puberulent with stem leaves numerous tubular florets, with several ray florets 
putative hybrid b ( H-b) long white puberulence with stem leaves numerous tubular florets, without ray floret 
putative hybrid ce (H-c) proximally glabrous with stem leaves numerous tubular florets, with several ray florets 


L. subspicata proximally glabrous with stem leaves numerous tubular florets, without ray floret 
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Table 2 Sites of variation and indels of ITS4-5 from related taxa 


Polymorphic sites 








Taxa j 1 157 1 2 254 2 3 4 4 4 5 5 6 
z 0 5 - 8 4 - 7 1 6 8 8 5 6 5 

0 159 0 9 255 4 9 4 0 5 7 4 6 

A1-5 A T G TCC C€ C CT G c c c G c c c 
A2-5 A c G TCC C C CT G c T T G c c c 
A3-1 A c G TCC C C cr A T T c G c c c 
B1-3 A c G TCC C c CT G c T T G c c c 
5-2 A c G TCC C€ C CT G T c c G c c c 
S1-S10, S12 A c G TCC C€ C CT G c T Y G c c c 
Al wo Y R ē YM Y - RY A c c c R Y Y c 
A2 w c R YM Y C RY G c Y Y R Y Y c 
A3 w c R YM Y - RY A c c c R Y Y Y 
Bl w c R YM Y C RÝ A c Y Y R Y Y Y 
c5 w c R YM Y - RY A Y c c R Y Y c 
V1-V4, V6-V10, VI2 T c A CU T - AC A c c c A T T T 
A1-1 T c A C T - AC A c c c A T T T 
A2-3 T c A CA T - AC A c c c A T T T 
A3-2 T c A CTA T - AC A c c c A T T T 
B1-6 T c A CTA T - AC A c c c A T T T 
C5-3 T c A CT T - AC A c c c A T T T 


Note; C+T=Y A+C=M A+T=W A+G=R G+T=K C+G=S 


Taxa: A, H-a; B, H-b; C, H-c; S, L. subspicata; D, L. vellerea. Numbers following taxon initials are sample numbers and clone numbers (if any). 


The strict consensus tree from parsimony analy- 
sis of 25 individuals was generated to show phyloge- 
netic relationship among them ( Fig. 2). The MP tree 
indicates that putative parents clustered restively 
branches with high bootstrap value. With cloned se- 
quences, each individual possessed one copy of one 
putative parent and clustered to parent branches ex- 
actly. 

2.3 Phylogenetic analysis based on sequences of 
cpDNA atpB-rbcL 

Direct sequencing of atpB-rbcL generated 740 
bp for L. vellerea and 739 bp for L. subspicata , exhib- 
ited 4 variation sites which included one indels ( Ta- 
ble 3). Among these sites, there were no polymor- 
phisms within L. subspicata individuals but one ex- 
ception in L. vellerea. All of putative hybrids had the 
same atpB-rbcL sequences with L. subspicata. 

The MP tree obtained from the cpDNA atpB- 
rbcL regions is shown in Fig. 3. The putative parental 
species formed two separate clusters. All the putative 


hybrids clustered together with L. subspicata. The 


other Ligularia species were clearly separated from 
the putative parental species and hybrids. Because 
chloroplast DNA is maternally inherited in Ligularia 
(Zhang et al., 2003), these putative hybrids had 
L. subspicata as the plastid donor. Hence unidirec- 
tional hybridization was detected between L. subspicata 
and L. vellerea. 


Table 3 Sites of variation and indels of chloroplast 


sequences from related taxa 


Polymorphic sites 








Taxa 4 5 6 
0 1 
Coa a 
S1-S10, S12 G T = c 
Al G T - c 
A2 G T - c 
A3 G T - c 
Bl G T - c 
c5 G T - c 
V1-V4, V6, V7, V9, 
V10, V12 j A 7 A 
V8 G A A A 


Taxa initials and numbering are as in Table 2 
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3 Discussion 

During our field research on Ligularia species 
more than decade, L. vellerea and L. subspicata are 
commonly found coexist, but no intermediate mor- 
phological individuals were discovered in these sym- 
patric areas except the one on Mt. Daxue. The puta- 
tive hybrids grow on an open space in a forest edge 
that was previously dominated by fir. Some trees 
were destroyed in road building and maintenance, 
and now grassland mainly grows the putative parental 
Ligularia species appeared on both sides of the road. 


This easily disturbed micro-habitat may have promo- 
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ted hybridization between these two species. Such 
special circumstances are also more likely to give 
rise to new species of hybrid origin ( Buerkle et al. , 
2000; Rieseberg et al., 2003; Gompert et al., 2006). 

Plant grow on wintry highland requires to evolve 
certain common adaptive mechanisms. Ligularia spe- 
cies often share the same pollinators in the highlands 
(Liu, 2002; Cao et al., 2008). When different spe- 
cies bloom at same place and time, nectar-collecting 
insects may transfer one species’ pollen to another spe- 
cies’ flowers casually. The flowering times of the sym- 
patric Ligularia species overlap, according to previous 


Ligularia platyglossa 
L. stenoglossa 

L. japonica 

L. atroviolacea 

L. thyrsoidea 

L. taquetii 

L. kongkalingensis 
L. curvisquama 

L. cyathiceps 

L. latihastata 

Vil 
v2 
V3 
V4 
V6 
v7 
V8 
v9 
V10 
V12 
Al-5 
A2-3 
A3-2 
B1-6 
C5-3 
Al-5 
A2-5 
A3-1 


L. vellerea 


Clones from 
putative hybrids 





B1-3 
C5-2 
S1 
S2 
S3 
S4 





S5 
S6 
S7 
S8 
S9 
S10 
S12 


L. subspicata 


Fig.2 Strict consensus tree from a parsimony analysis of ITS4-5 sequences from putative hybrid individuals and their putative parental species. 


Taxa: V, L. vellerea; S, L. subspicata; A, putative hybrids a; B, putative hybrids b; C, putative hybrids c. Numbers following 


taxon initials are sample numbers and clone numbers (if any). Bootstrap percentages >50% are shown above the branches 
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S4 
S6 


L. subspicata 
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S8 
S9 
S10 
S12 
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A2 








A3 putative hybrids 


B1 
C5 
Vi 
V2 
V3 
V4 
V6 
V7 
v9 
V10 
V12 
V8 


L. vellerea 














L. franchetiana 





L. oligonema 








L. anoleuca 
L. veitchiana 
L. kanaitzensis var. subnudicaulis 


Fig.3 Maximum parsimony tree inferred from atpB-rbcL sequences to show the relationships between the putative hybrids and putative parents. 


Bootstrap percentages >50% are shown above the branches. Taxon initials and numbering are as in Fig. 2 


researchers and our field observation for years. They 
flower in late June to August, but their flowering 
times differ somewhat; the anthesis L. vellerea is a little 
earlier than L. subspicata. 

Our ITS4-5 data definitely support the hypothesis 
that the morphologically intermediate individuals are 
the progeny of hybridization between L. vellerea and 
L. subspicata. According to the results of cpDNA analy- 
sis, we can also conclude that the hybridization is uni- 
directional in the form of L. subspicata Ẹ X L. vellerea 3. 

Three types of putative hybrids were assigned 
according to their different morphology. We thought 
this hybridization might bidirectional and parental 
species might more than the putative two. All of hy- 
brids we studied were more close to parental species 
of L. subspicata. But with some different morphologi- 
cal characters in detail, the first hybrids type H-a 
had L. vellerea florets component part with several 


ray florets except tubular florets in one capitulum, 


H-b had hairs on stem and H-ce had ray florets. 
Though they had same plastid donor, they might in- 
herit different copies of nuclear genetic material 
which leads to different phenotype. Another possibil- 
ity is F,, Fns or backcross appeared and hybrid 
swarm formed at this population. Because materials 
were limited in the first collection, new and long 


term observation was needed for further studying. 
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Appendix; List of accessions of Ligularia included in the phylogenetic analysis of nuclear rDNA ITS and atpB-rbcL sequences. 


GenBank numbers for all accessions are provided. 


Accessions , GenBank No. 


ITS: Ligularia platyglossa, AY458841; L. stenoglossa, AB523365; L. japonica, AY458825; L. atroviolacea, DQ272322; 
L. thyrsoidea, DQ272320; L. taquetii, GU323282; L. kongkalingensis, DQ272343; L. curvisquama, DQ272324; L. cyathiceps , 


DQ272328; L. latihastata, DQ272334. 


atpB-rbcL; L. franchetiana, AB375447; L. oligonema, AB523847; L. anoleuca, AB557885; L. veitchiana, AB557887; 


Ligularia kanaitzensis var. subnudicaulis, AB281495. 


